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Abstract — Three-dimensional Particle-in-Cell (PIC) simu- 
lations with the code QuickPIC are used to illustrate the 
typical accelerating structures associated with the interac- 
tion of an intense laser beam with an underdense plasma in 
the blowout regime. Our simulations are performed with an 
externally injected electron beam, positioned in the region 
of maximum accelerating gradients. As the laser propagates 
in the plasma, almost complete electron cavitation occurs, 
leading to the generation of accelerating fields in excess of 
1 GeV/cm. 
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The recent developments on infrared (Ao ~ A«n), ultra- 
intense (> 10 19 W/cm ), ultra-short (< 30 fs), and tightly 
focused (~ 20 /jm) laser beams have opened the possibil- 
ity to explore new scenarios on laser-plasma interactions, 
namely electron acceleration in the bubble [l] or blowout [2] 
regime of the laser wakefield accelerator (LWFA) |3j. In 
many LWFA experiments [4], [5], [6], [7] the laser intensity 
is sufficiently high for the radiation pressure of the laser 
to cause complete electron cavitation. Since ions are es- 
sentially stationary during a plasma wave period, a strong 
space charge is formed that pulls electrons back to the axis. 
The wake then assumes a spherical shape [l] , [2] , with lin- 
ear accelerating and focusing forces. In addition to the 
non-linear plasma response, the laser self-modulations also 
play a key role in LWFA experiments. Thus, the overall 
laser-plasma dynamics is strongly coupled, and a complete 
theoretical understanding of typical LWFA experiments is 
beyond the capabilities of current pure theoretical models. 

Simulations play an important role in this matter, by giv- 
ing deeper insights on the physics of the LWFA, and being 
also valuable in the modeling and design of new experi- 
ments. In the pursuit of higher energy gains, lower plasma 
densities are required. Although the peak accelerating field 
is reduced for lower plasma densities with -E a ccei oc ■ s pn e , 
where n e is the plasma density, the dephasing/pump- 
depletion lengths increase with L acco i cx n e 3 ^ 2 . There- 
fore, the maximum energy gain increases with W mm ~ 
-Eaccei-kaccei oc n~ x . The use of lower densities, however, 
constitutes a major difficulty to LWFA simulations using 
standard full PIC codes: it requires both finer grids to 
correctly resolve the laser wavelength, and more simula- 
tion time steps to model the increased interaction length. 
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Therefore, the use of reduced PIC models to examine the 
laser-plasma interaction in the LWFA is becoming increas- 
ingly important. In this paper, we use the reduced PIC 
code QuickPIC (8| to illustrate the properties of the wake 
created by an intense laser pulse in the presence of an ex- 
ternal electron beam. QuickPIC works under the quasi- 
static approximation (QSA) [9]. The QSA assumes that 
the typical time for the laser evolution is much larger than 
the typical time for the plasma response, precluding, how- 
ever, the physics associated with self-injection. This as- 
sumption is very well verified in the LWFA. In QuickPIC, 
the plasma response is then calculated according a quasi- 
static field solver for each transverse plasma slice and the 
laser is advanced through the ponderomotive guiding cen- 
ter approximation. With this approximation larger time 
steps can be used for the laser evolution, thus leading to 
computational time savings of more than three orders of 
magnitude in comparison to standard full PIC codes. 

The three-dimensional structures associated with the 
laser propagation in the blowout regime are illustrated in 
Figure [l] We use a laser pulse with a central wavelength 
Ao = 0.8 //m, focused to 20 /jm, with a duration of 30 fs, 
and a peak normalized vector potential ao = 4. The laser 
pulse propagates in a plasma with density n e — 10 18 cm~ 3 . 
A bi-gaussian electron beam is externally injected, with a 
transverse width of 2 fim, duration of 20 fs, charge of 1 nC, 
and traveling with a relativistic factor 7 = 10 4 . The sim- 
ulation box moves at the speed of light and it is resolved 
with 256 x 256 x 512 cells for the transverse and longitudinal 
directions respectively, with 16 particles per cell. 

The plasma electrons are almost completely evacuated 
from the region where the laser is sitting (cf. Fig. [l]-(a)). 
As they return to the axis (cf. Fig. [T]-(b)), the plasma 
electrons feel the repulsive radial Coulomb force associated 
with the externally injected beam, leading to the flatten- 
ing of the linear focusing force (cf. Fig. [T]-(c)). As a con- 
sequence, the wake does not fully close at the end of the 
first plasma period. Similarly, in the longitudinal direction 
the accelerating gradient is also flattened by the presence 
of the electron beam (cf. Fig. [T]-(d)). 

In conclusion, we have presented figures of three- 
dimensional waves in the blowout regime. These im- 
ages clearly illustrate the fundamental structures associ- 
ated with the LWFA, and allow a better understanding of 
the physics associated with these experiments. 
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Fig. 1 

(a) Isosurfaces of the plasma and externally injected beam charge density (green). The isosurfaces of the laser fields are 
shown in the front of the box (gold). the laser propagates in xl. the projection corresponds to the plasma and beam 

DENSITIES, AND TO THE LASER INTENSITY, (b) INTENSITY ISOSURFACES AND FIELDLINES (HIGHER CURRENT VALUES ARE IN RED AND LOWER 
CURRENT VALUES ARE IN BLUE) OF THE PLASMA CURRENT IN THE BLOWOUT REGION. THE PROJECTION CORRESPONDS TO THE INTENSITY OF 

THE PLASMA CURRENT, (c) ISOSURFACES OF THE FOCUSING FORCE, WITH THE CORRESPONDING PROJECTIONS, (d) ISOSURFACES OF THE 
LONGITUDINAL ACCELERATING GRADIENT. THE PROJECTIONS ALSO CORRESPOND TO THE ACCELERATING WAKE FIELD (HIGHER FIELD VALUES 

ARE IN RED AND LOWER FIELD VALUES ARE IN BLUE) 
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